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Abstract. We present the results of a mid-infrared 
(A ~ 7 fixsx) imaging survey of a sample of 24 starless dense 
cores carried out at an angular resolution of 6" with the 
ISOCAM camera aboard the Infrared Space Observatory 
(ISO). The targeted cores are believed to be pre-stellar 
in nature and to represent the initial conditions of low- 
mass, isolated star formation. In previous submillimeter 
dust continuum studies of such pre-stellar cores, it was 
found that the derived column density profiles did not fol- 
low a single power-law such as Nn 3 k F 1 throughout 
their full extent but flattened out near their center. These 
submillimeter observations however could not constrain 
the density profiles at radii greater than ~ 10000 AU. 
The present absorption study uses ISOCAM's sensitiv- 
ity to map these pre-stellar cores in absorption against 
the diffuse mid-infrared background. The goal was to de- 
termine their structure at radii that extend beyond the 
limits of sensitivity of the submillimeter continuum maps 
and at twice as good an angular resolution. Among the 24 
cores observed in our survey, a majority of them show deep 
absorption features. The starless cores studied here all 
show a column density profile that flattens in the center, 
which confirms the submillimeter emission results. More- 
over, beyond a radius of ~ 5000 — 10000 AU, the typ- 
ical column density profile steepens with distance from 
core center and gets steeper than Nu 2 oc f , until it 
eventually merges with the low-density ambient molecu- 
lar cloud. At least three of the cores present sharp edges at 
R ~ 15000 — 30000 AU and appear to be decoupled from 
their parent clouds, providing finite reservoirs of mass for 
subsequent star formation. 
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Although our general understanding of low-mass star for- 
mation has made significant theoretical and observational 
progress in the past three decades (see, e.g., the volumes 
by Levy & Lunine 1993 and Mannings, Boss, & Russell 
2000 for comprehensive reviews) , the earliest stages of the 
star formation process remain poorly known. These stages 
immediately preceding and immediately following the on- 
set of protostellar collapse are of great interest since to 
some extent they must govern the origin of stellar masses, 
i.e., the stellar initial mass function. 

According to present ideas on isolated star forma- 
tion (e.g. Shu, Adams, & Lizano 1987), the first evolu- 
tionary phase in the path from molecular cloud to main 
sequence star involves the formation of gravitationally- 
bound, starless dense cores (e.g. Myers 1999). These "pre- 
protostellar" (or "pre-stellar" for brevity) cores (cf. Ward- 
Thompson et al. 1994) are thought to be initially sup- 
ported against their self-gravity by magnetic and/or tur- 
bulent pressure, and to progressively evolve toward higher 
degrees of central concentration through ambipolar dif- 
fusion (e.g. Mouschovias 1991) and/or the dissipation of 
turbulence (e.g. Nakano 1998). At some not yet well un- 
derstood point, the cores become unstable and collapse 
dynamically to form accreting protostars which quickly 
become detectable in the radio range as Class objects 
(Andre, Ward-Thompson, & Barsony 1993) and subse- 
quently in the infrared as Class I sources (eg. Lada 1987). 
Additional complications arise in regions of multiple star 
formation (such as the p Ophiuchi central cloud) where 
external triggers rather than self-initiated ambipolar diffu- 
sion may be responsible for cloud fragmentation and core 
formation (e.g. Loren & Wootten 1986; Boss 1995; Whit- 
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worth et al. 1996; Motte, Andre, & Neri 1998 - hereafter 
MAN98). 

A good knowledge of the initial conditions for rapid 
protostellar collapse in star-forming cores is crucial to get 
at a proper theoretical description of early protostar evo- 
lution. In particular, recent studies show that the mass- 
infall rate during the protostellar phase is quite sensitive 
to the core radial density profile at the onset of dynamical 
collapse. When the initial density profile is flatter than 
p oc r -2 near core center, then the mass-infall rate is ex- 
pected to feature a strong initial peak, followed by a sharp 
decline, at the beginning of the main accretion phase (e.g. 
Foster & Chevalier 1993; Henriksen, Andre, & Bontemps 
1997). The outer core density profile also plays a funda- 
mental role by determining whether the mass reservoir 
ultimately available for star formation is effectively finite 
or infinite. If the radial density profile approaches p oc r~ 2 
up to large radii as in the standard paradigm of isolated 
star formation (Shu et al. 1987), then the mass reservoir is 
much larger than a typical stellar mass and feedback pro- 
cesses such as outflow/inflow interactions (e.g. Velusamy 
& Langer 1998) must eventually stop accretion onto the 
central protostar. On the other hand, if the initial core 
density profile becomes steeper than p oc r~ 3 beyond some 
finite radius, then the mass reservoir is limited and the 
mass accretion rate will naturally drop to negligible val- 
ues after some time. In this case, the final stellar mass 
may be essentially determined at the pre-collapse stage. 

Direct observational constraints on the radial density 
structure at the onset of protostellar collapse have been 
obtained through (sub)millimeter dust continuum map- 
ping of a large number of starless cores and Class O/Class 
I protostellar envelopes with the JCMT and IRAM 30 m 
telescopes (e.g. Ward-Thompson et al. 1994; Andre, Ward- 
Thompson, & Motte 1996 - hereafter AWM96; Ward- 
Thompson, Motte, & Andre 1999 - hereafter WMA99). 
General trends have been found: while protostellar en- 
velopes are always highly centrally condensed, with ra- 
dial density profiles consistent with p oc r -1 ' 5 or p oc r~ 2 
(e.g. Ladd et al. 1991, Motte & Andre 1999), starless cores 
have radial profiles that flatten out near their centers and 
become much flatter than p oc r~ 2 at radii less than a 
few thousand AUs (see Andre, Ward-Thompson, & Bar- 
sony 2000 for a review) . It is however necessary to bear in 
mind the limitations of (sub) millimeter continuum emis- 
sion maps: since the dust emission from the outer parts of 
cores is intrinsically very weak, (sub)millimeter mapping 
is mostly sensitive to the inner density structure of star- 
less cores. The outer density structure of cores and the 
possible presence of edges arc difficult to constrain from 
(sub)millimcter observations. Fortunately, starless cores 
can also be studied in absorption at mid-IR wavelengths 
(e.g. Abergel et al. 1996, Egan et al. 1998), and this tech- 
nique turns out to be more sensitive to the outer parts of 
cores. Another advantage of the mid-IR absorption tech- 
nique over (sub)millimeter emission mapping is that the 



column density structure of the absorbing material may 
be derived without any assumption about the dust tem- 
perature distribution. 

In order to gain further, independent insight into cloud 
core structure prior to protostar formation, we undertook 
dedicated mid-IR imaging observations toward a broad 
sample of starless dense cores with the ISOCAM mid-IR 
camera aboard the Infrared Space Observatory (ISO). The 
present paper discusses the results of this mid-IR imaging 
program, whose immediate goal was to detect the sample 
cores in absorption against the diffuse mid-IR background 
produced by the mean galactic interstellar radiation field 
(e.g. Mathis, Mezger & Panagia 1983) and/or arising from 
the envelopes of the parent molecular clouds (e.g. Bernard 
et al. 1993). 

The layout of the paper is as follows. In Sect. 2, we give 
the properties of the core sample and explain observational 
details. In Sect. 3, we show our ISOCAM absorption maps 
and analyze the associated absorption profiles in terms of 
column density structure. In Sect. 4, we compare our re- 
sults with various theoretical models of core structure and 
discuss possible implications for our understanding of the 
pre-stellar stage of star formation. Our main conclusions 
are summarized in Sect. 5. 

2. Observations and Data Reduction 

2.1. Starless core sample 

The 24 starless fields of the present survey were selected 
from various dark cloud and nebula catalogs (e.g. Parker 
1988, Benson & Myers 1989, Bourke, Hyland & Robin- 
son 1995, Schneider & Elmegreen 1979), on the grounds 
that they contained no IRAS point sources associated 
with a young stellar object and had an estimated mid- 
infrared (mid-IR) background stronger than ~ 1 MJy/sr 
at ~ 7 pm. We also included a map of the dark cloud 
Oph D (also called L1696A) situated in the p Ophiuchi 
complex, and already observed in absorption by ISOCAM 
(Abergel et al. 1998). The background at 7 pm was orig- 
inally estimated from IRAS images, as were the 7 /xm 
flux densities of the IRAS sources appearing in the fields. 
Fields containing IRAS point sources with estimated flux 
densities larger than 1.5 Jy at 7.75 pm (LW6 filter [7- 
8.5 /im]) or larger than 0.7 Jy at 6.75 pm (LW2 filter [5- 
8.5 pm]) had to be shifted away from the sources and/or 
the integration time had to be shortened in order to avoid 
saturating the ISOCAM detectors (see 2.2 below). The 
selected regions span a wide range of properties (e.g. mor- 
phologies, environments), but are all located in nearby 
molecular cloud complexes (distance < 500 pc) either ac- 
tively forming stars or supposed to be future hosts of star- 
formation activity. 

The properties of the observed cloud cores are summa- 
rized in Table 0. Column 1 lists the names of the sources 
as found in the dark cloud catalogues, columns 2 and 3 
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Table 1. List of starless cores observed with ISOCAM 



Source R. A. (2000) Dec. (2000) Dist. Filter Map RMS <Imir > ll odi Absorp. References 









(PC) 




(MJy/sr) 


(MJy/sr) 


(MJy/sr) 


Contrast 

/(V \ 

(%) 




L1517B* 


04 ,l 55 m 05':2 


30°38'39" 


140 


bW2 


0.07 


4.4 


3.8 


5 


(4),(6),(13) 


L1512* 


05 h 04 m 09!7 


32°43'09" 


140 


LW6 


0.15 


12.4 


- 


3 


(4),(5),(6),(H),(15) 


L1544* 


05 h 04 m 18 s 1 


25°11'08" 


140 


LW2 


0.05 


4.8 


4.2 


8 


(4),(11),(14),(15),(16) 


L1582A 


05 h 32 m 03!0 


12°30'16" 


400 


LW2 


0.16 


8.6 


3.7 


44 


(4),(6) 


L1672 


05 h 54 m 24 s 


01°58'53" 


200 


LW2 


0.19 


4.7 


3.4 


31 


(10) 


BHR78 


12 fe 36 m 18!7 


-63°12'37" 


200 


LW6 


0.61 


17.2 


- 


11 


(3) 


BHR111 


15 h 42 m 19 s 7 


-52°47'54" 


250 


LW6 


0.41 


19.5 


10.0 


6 


(3) 


R7** 


16 h 21 m 4S!9 


-23°42'51" 


160 


LW6 


0.92 


- 


- 


- 


(8) 


OphD** 


16 h 28 m 3ff4 


-24° 18' 29" 


160 


LW2+ 


0.3 


14.8 


5.7 


28 


(1),(8),(9) 


R53" 


16 fe 31 m 44!l 


-24° 48' 12" 


160 


bW2 


0.2 


11.4 


- 


20 


(8) 


SA187 




-44°53'40" 


200 


LW6 


0.92 


- 


- 


6 


(12) 


bl709A** 


16' l 32 m 42=l 


-23°54'09" 


160 


LW2 


0.17 


12.8 


10.5 


17 


(4),(6),(8),(11) 


R60" 


16 fc 33 m lK8 


-24°41'30" 


160 


LW2 


0.13 


10.1 


- 


7 


(8) 


L1709C" 


16 h 33™ 53" 4 


-23° 42' 32" 


160 


LW2 


0.10 


8.2 




16 


(4),(6),(8),(11) 


R63" 


16 fc 33 m 58S3 


-24°30'03" 


160 


LW2 


0.14 


9.3 


5.0 


16 


(8) 


L1689B" 


16 h 34 m 4(n 


-24°37'00" 


160 


LW2 


0.05 


8.3 


4.3 


19 


(2),(4),(8),(H),(15) 


L204B 


16 h 47 ™ 32 s 


-11°59'15" 


170 


LW2 


0.11 


15.7 




4 


(4),(H) 


B68 


17 h 22 m 34 s 7 


-23°47'50" 


200 


LW2 


0.19 


11.5 


8.8 


2 


(4),(10) 


BHR140 


17 ft 22 m 53l9 


-43°22'13" 


400 


LW6 


0.34 


10.2 


6.8 


9 


(3) 


b310 


18 h 07 m u s 9 


-18°21'35" 


200 


LW6 


0.54 


30.0 


7.1 


39 


(11) 


L328 


18 ft 17 m 0(K4 


-18°01'52" 


200 


LW6 


0.7 


42.9 


7.4 


40 


(11) 


L429 


18H7 m 06M 


-08°15'51" 


200 


LW6 


0.2 


8.6 


6.8 


17 


(11) 


GF5 


18 h 39 m 16!4 


-06°38'15" 


200 


LW6 


2.0 


75.2 


7.2 


30 


(13) 


B133 


19 ft 09 m 08l8 


-06°53'20" 


400 


bW6 


0.13 


7.9 




3 


(4),(7) 



(l)=Abergel et al. 1996; (2)=AWM96; (3)=Bourke, Hyland & Robinson 1995; (4)=Benson & Myers 1989; 
(5)=Caselli, Myers, Thaddeus 1995; (6)=Hilton & Lahulla 1995; (7)=Keene 1981; (8)=Loren 1989; (9)=Loren et 
al. 1990; (10)=Leung, Kutner & Mead 1982; (ll)=Parker 1988; (12)=Sandqvist & Lindroos 1976, Sandqvist 
1977; (13)=Schneider & Elmegreen 1979; (14)=Tafalla et al. 1998; (15)= Ward- Thompson et al. 1994; 
(16)=Williams et al. 1999 

^Oph D was also mapped in the LW3, LW6 and LW9 filters (see text and Fig. ^cd). The zodiacal intensity is 45 
MJy/sr in LW3. 

^Values of l zo di come from the zodiacal light model of Kelsall et al. (1998). Typical uncertainty is ~ 20%. 
* Cores belonging to Taurus molecular complex. 
"Cores belonging to p Ophiuchi molecular complex. 



list the J2000 coordinates of the centre of the ISOCAM 
field. The image centres are sometimes shifted from their 
original positions as given in the dark cloud catalogues 
(see above). Column 4 gives the distances to the clouds. 
For the clouds not belonging to the p Ophiuchi (adopted 
distance d = 160 pcQ) and Taurus (d — 140 pc) molecular 
complexes, we used the values from the study of Hilton 

6 Lahulla (1995) or the results of Dame et al. (1987) 
on the locations of the major (giant) molecular clouds in 
the Galaxy to give an estimate of the distances by relat- 
ing the cores to known molecular complexes according to 
their galactic coordinates. Column 5 gives the ISOCAM 
filter(s) we used for the observations of each core. Column 

7 lists the mean mid-IR intensity measured off the absorp- 
tion features. Included in the values listed in column 7 is 



the zodiacal light emission (given in Column 8). Column 
9 lists the relative contrast between the intensity of the 
main absorption feature and the mean mid-IR intensity 
(column 7) , in percentage of the mean mid-IR intensity q 
Finally, column 10 gives the references of the dark 
cloud catalogues, from which our fields were selected, as 
well as the references of a few studies containing interest- 
ing information on the sources. From Table ^, it can be 
seen that the majority of the fields we observed contain ab- 
sorption features. Figure ^ shows eleven cores with deep 
absorption (i.e., L1544, L1582A, Oph D, R53, L1709A, 
R63, L1689B, L310, L328, L429, GF5 - see Section O)). 



1 although recent studies suggest d ~ 120—140 pc 
Knude & Hog (1998) and de Zeeuw et al. (1999). 



2 Since the intensity in column 7 includes the contribution 
of the zodiacal light emission, the contrast in column 9 is not 
truly representative of the actual core absorption strength. 
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2.2. ISOCAM observations 

2.2.1. Observational parameters 

The dense cores were mapped between April 1997 and 
March 1998 using the LW2 (5-8.5 pcm) and LW6 (7-8.5 
(im) filters of ISOCAM (Cesarsky et al. 1996), the spectro- 
imager aboard ISO (Kessler et al. 1996), centered at 6.75 
/im and 7.75 /im respectively (except for Oph D, see be- 
low). These wavelengths approximately correspond to a 
minimum in the dust extinction curve (see the model of 
Draine & Lee 1984 for diffuse clouds), providing a good 
opportunity to peer deep inside dark clouds. We used a 
pixel field of view of 6", which also roughly corresponds 
to the FWHM size of the ISOCAM point spread func- 
tion (PSF) at ~ 7 /im. For a "standard" observation, an 
area of 10.7' x 11' was covered for each core with a raster 
map of 6 x 4 images and an overlap of 17 pixels in right 
ascension and of 6 pixels in declination. At each raster 
position, about 18 single read-outs were taken with an in- 
tegration time per read-out of 2.1s. About 50 stabilisation 
read-outs had to be added at the beginning of each obser- 
vation to allow the detectors to lose memory of the previ- 
ous observation. For a few dense cores (L1672, SA187, and 
GF5) , the field contained strong IRAS point sources and a 
shorter (0.28s) integration time per read-out was adopted 
in order to avoid saturation on the camera detectors. In 
these cases, the number of read-outs was 44 at each raster 
position and the total area mapped was 19.7' x 18.8' for 
L1672 and SA187. In addition to the existing LW2 (5-8.5 
H and LW3 (12-18 /tin) maps of Abergel et al. (1996, 
1998) with 6" pixels, we took two new images of Oph D 
in the filters LW6 (7-8.5 /im) and LW9 (14-16 /im) with 
a pixel field of view of 3". The Oph D field was covered 
by a raster of 6 x 3 images with an overlap of 22 pixels 
in both directions. The raster was tilted by 20° with re- 
spect to the Celestial North in order to follow the core's 
elongation (see Fig. 2d of MAN98 and Fig. 1 of Abergel et 
al. 1998). For these images of Oph D, 15 exposures were 
taken at each raster position with an integration time of 
10 sec per read-out. The total area covered was 2.6' x 4.1'. 

2.2.2. Data reduction 

The data were reduced using the CIA|^| (CAM Interactive 
Analysis) software package. First, glitches due to cosmic 
rays were removed using a spatial and temporal filtering 
algorithm (Starck et al. 1999), then a dark current im- 
age was subtracted from all the frames. Because of the 
time lag in the response of the detectors, transient effects 
had to be corrected for. For this purpose, we used the al- 
gorithm developed by Abergel et al. (1996). Finally, the 
images were divided by a flatfield image to correct for the 



3 'CIA' is a joint development by the ESA Astrophysics Divi- 
sion and the ISOCAM Consortium. The ISOCAM Consortium 
is led by the ISOCAM PI, C. Cesarsky. 



non-uniform response of the detectors over the array. Most 
images were made up of over 500 frames that covered a 
large surface area. This enabled us to use a median image 
of all the frames as a flatfield image. Flatfielding is here 
a critical operation as we are looking for intensity varia- 
tions in dark objects and any flux mismatch between two 
contiguous raster images appears as a step in an inten- 
sity profile. Additional pixels had to be masked where the 
trails of glitches still affected the response of the detectors 
or induced a downward step in the flux level. 

2.2.3. Noise estimation 

In order to estimate the rms noise levels in our final im- 
ages, we first derived the rms level of the temporal fluc- 
tuations for each pixel of a given raster map, taking into 
account the overlap of the individual raster frames (each 
sky position in the map was observed on average in two in- 
dependent individual frames). The resulting average rms 
noise per pixel depends on the level of the diffuse mid- 
IR emission in the raster, and thus varies from core to 
core. Typical values are ~ 0.05 MJy/sr for regions with 
low mid-IR background (e.g. Taurus, Eastern part of the 
p Oph complex) and ~ 0.1 MJy/sr for regions with high 
mid-IR background (e.g. L328 and L310). This temporal 
rms noise however does not take flatfield noise (estimated 
to be ~ 5%) into account and could be underestimating 
the total amount of noise. Another approach to estimating 
the noise in an image is by calculating the rms noise level 
of the spatial fluctuations over regions showing very little 
structure. By doing so, we a priori overestimate the noise 
level since we do not discriminate between true instrumen- 
tal noise in the image and small-scale spatial structure in 
the (emitting and/or absorbing) source. The rms values 
listed in Ta ble [l] and the error bars used in the profiles 
of Sect. 3.4 below correspond to such spatial estimates 
and should thus be taken as upper limits. Note that in 
the densest areas, we can distinguish between small-scale 
structure and instrumental noise or mid-IR emission fluc- 
tuations by checking for the presence on small scales of 
cold, IR-absorbing condensations in high-resolution mil- 
limeter continuum maps (see 2.3 below). 



2.3. Millimeter observations 

In order to check that the dark patches seen in the ISO- 
CAM images (see Fig. Q) were actually due to the pres- 
ence of cold absorbing dense cores/condensations, and 
not to small-scale fluctuations of the mid-IR background 
and/or foreground, we carried out 1.3 mm continuum ob- 
servations of a selected subset of our fields (see Table @). 
We used the IRAM 30 m telescope equipped with the 
37-channel MPIfR bolometer array in March 1998 and 
March/April 1999 to map the 1.3 mm dust emission of 
those fields that presented clear and compact absorption- 
like features in the ISOCAM images. On-the-fiy 1.3 mm 
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maps were taken in the dual-beam raster mode with a 
scanning velocity of 4"/sec, a spatial sampling of 2" in 
azimuth and 4" in elevation, and a chopping frequency of 
2 Hz. The chop throw was set to 60" for maps larger than 
300" and 45" for smaller maps. The (FWHM) beam size of 
the telescope was measured to be approximately 11"-12" 
on maps of Uranus. The zenith atmospheric optical depth, 
monitored by 'skydips' every 1-2 hours, was between ~ 
0.15 and ~ 0.5. The final co-added maps are about 7' x 7' 
in size and consist of ~ 5-6 individual coverages taken at 
different hour angles to improve the reconstruction tech- 
nique. The data were reduced with the IRAM software 
for bolometer arrays ("NIC"; cf. Broguiere, Neri & Sievers 
1995), which uses the EKH restoration algorithm (Emer- 
son, Klein, & Haslam 1979). The typical rms noise was 
~ 8 mJy/ll"beam. 

The same dense cores listed in Table || were also ob- 
served in the C 18 O(l-0) molecular line at the IRAM 30 m 
telescope, in order to obtain estimates of the column den- 
sity in the outer parts of the cores. We used the position 
switching mode along cuts of about 10' length sampled ev- 
ery 20" through each of the cores, and oriented along and 
perpendicular to the axis of the core. The OFF position, 
which was checked to be free of emission, was typically 
taken between 15' and 30' away from the nominal source 
position. Our goal was to get good signal-to-noise ratio in 
the outer parts of the core in order to estimate the H2 
column density, Njj 2 , as far from the core centre as possi- 
ble (typically ~ 200"). The total (ON+OFF) integration 
time was 2 minutes on average and up to 10 or 25 minutes 
in the outer parts of the cores in order to reach an rms 
sensitivity better than ~ 0.1 K km s _1 . The spectral res- 
olution was 20 kHz and the observing frequency 109 GHz. 
Additional N 2 H+(l-0), H 2 CO(2-l), and DCO+(2-l) spec- 
tra were taken at a few positions for selected cores. 

During these continuum and line runs at the 30 m 
telescope, the pointing, as well as the focus, was checked 
every ~ 1 hour and found to be better than ~ 5" in both 
azimuth and elevation. 

3. Results and analysis 

3.1. ISOCAM images 

Out of the 24 dark cloud fields in our programme, 23 show 
absorption-like features in the mid-IR. 

The fields presenting such features can be divided into 
two broad categories, those that show deep ( £ 10% - 
cf. column 9 of Table [l]), compact absorption features 
likely associated with dense cores, and those presenting 
weaker ( ^ 10%) and more diffuse absorption whose na- 
ture is also more uncertain. Figure [j] presents the ISO- 
CAM images of 11 fields with strong absorption features, 
to which we added the widely studied core L1544 (e.g. 
Ward-Thompson et al. 1994, Tafalla et al. 1998, Williams 
et al. 1999, Ohashi et al. 1999), and from which we with- 
drew L1709C (the absorption features lying in the corners 



of the image, making the analysis difficult). The remain- 
ing fields are presented in Fig. ||. Figure || shows examples 
of cuts through two of the deep absorption features (for 
L328 and L1689B). The morphologies encountered in the 
fields with compact absorption vary from spherical cores 
like L328 to more elongated, filamentary structures (e.g. 
L1689B, GF5). 

3.2. Comparison with the millimeter continuum maps 

All of the cores observed at 1.3 mm (Table |) present ex- 
tended, weak dust continuum emission which we were able 
to map. For all cores, the peak of the mid-IR absorption 
does correspond with a peak in the millimeter emission. 
This is illustrated in Fig. |^ which shows the mid-IR ISO- 
CAM absorption contours superimposed on the 1.3 mm 
continuum maps of the eight cores positively detected at 
1.3 mm (the maps of Oph D and L1544 are also shown in 
MAN98 and WMA99, respectively). 

The good agreement between the respective images of 
these cores in both types of data confirms that the features 
we see with ISOCAM are due to absorption by the (rather 
large) column density of cold dust traced at 1.3 mm, and 
not by fluctuations in the mid-IR background. Further- 
more, the similar morphologies seen in both wavebands 
suggest that the mid-IR absorption images trace the same 
dust as the millimeter emission maps, which justifies a de- 
tailed comparison of the column density profiles derived 
from both tracers. 

3.3. Modelling of the mid-infrared absorption 

In order to derive information about the column density 
structure from the mid-IR intensity, we have used a simple 
geometrical model: each dense core is embedded within a 
lower-density parent molecular cloud. On the line of sight 
to the core, we measure i) a background intensity arising 
from the rear side of the parent cloud that is attenuated by 
the absorption from the core, and ii) a foreground intensity 
arising from the front side of the parent cloud. We can thus 
express the mid-IR intensity along the line of sight as: 

1(f) = [< hack > +SI ba ck(f)} ■ e- T ^ 

+ [< I fore > +SIfore(r)], (1) 

where hack and I fore ar e the background and foreground 
intensities respectively, and T\(r) is the dust opacity at a 
projected radius f from core center and at observed wave- 
length A. The foreground emission, If ore (f), includes the 
contribution from the zodiacal light emission, I zo di- 

The mid- infrared background emission, Ibacki^), is 
thought to arise from very small grains presenting aro- 
matic emission features (Leger & Puget 1984, Allaman- 
dola, Tielens & Barker 1989, Boulanger et al. 1996). These 
grains/molecules are excited by the interstellar far-UV 
(FUV) radiation field and re-emit this energy as a group of 
spectral lines/bands in the mid- infrared, the Unidentified 
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Fig. 1. ISOCAM images of the sample of 11 dense cores presenting strong absorption features (see text). The cores 
appear as dark patches against a brighter background emission. The ISOCAM psf is ~ 6" (FWHM) at ~ 7 /xm; 
the pixel size is also 6" (3" in the case of Oph D). W hite c rosses and triangles show the core centers and positions 
where the C 18 column density was measured (section 3.3.2| ), respectively. The dotted lines/sectors illustrate the way 
averages were performed (section [^4|). The greyscale is indicated in MJy/sr at the right of each image. Contour levels 
are: a) for L1544, 4.55, 4.65, 4.75 MJy/sr; b) for L1582A, 5, 5.5, 6 MJy/sr; c) for Oph D LW6 (7.75 /mi), 21, 22, 23, 
24 MJy/sr; d) for Oph D LW9 (15 (mi), 74, 75, 76 MJy/sr; e) for R53, 9.5, 9.75, 10 MJy/sr; f) for L1709A, 24, 25, 26 
MJy/sr. 
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Fig. 1. (cont'd) Contour levels: g) for R63, 7.75, 8, 8.25 MJy/sr; h) for L1689B, 7, 7.5, 8 MJy/sr; i) for L310, 20, 22.5, 
25, 27.5 MJy/sr; j) for L328, 30, 34, 38 MJy/sr; k) for L429, 7.3, 7.7, 8.1, 8.5 MJy/sr; 1) for GF5, 57, 62, 67 MJy/sr. 
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Fig. 2. ISOCAM images of the remaining cores. The fainter absorption features led us to enhance the contrast of the 
images, which resulted in poorer image appearance (the flatfield noise is clearly visible in some cases). Contour levels 
are: a) for L1517B, 4.3 MJy/sr; b) for L1512, 12.2 MJy/sr; c) for L1672, 3.5, 4 MJy/sr; d) for BHR78, 15, 16, 17; e) 
for BHR111, 17, 18, 19 MJy/sr; f) for R7, 20, 25, 30, 35 MJy/sr. 
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Fig. 2. (cont'd) Contour levels: g) for R60, 9.1, 9.3 MJy/sr; h) for L1709C, 6.9, 7.4 MJy/sr; i) for L204B, 12.2 MJy/sr; 
j) for B68, 5.4, 5, 6 MJy/sr; k) for BHR140, 9, 9.5 MJy/sr; 1) for B133, 7.4 MJy/sr. 
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Fig. 3. (a) North-South intensity cut through the A = 7.75 fmx ISOCAM image of L328 (see Fig. fD), at a ~ 18' l 17 m 0ffi 
(b) North-South intensity cut through the A = 6.75 fj,ra ISOCAM image of L1689B at a ~ 16%4 m 50! Both cloud 
cores clearly show up in absorption against the diffuse mid-infrared background. 



Infrared Bands (UIBs). The ISOCAM LW2=[5.0, 8.5 H 
and LW6=[7.0, 8.5 fxm] filters that we used for our obser- 
vations include the 6.2 /jm and 7.7 /im emission features 
of these UIBs. 

According to recent studies, the small particles respon- 
sible for the UIBs are distributed throughout molecular 
clouds but because of high interstellar extinction in the 
FUV band, they are ususally only excited in the outer re- 
gions where A v ^ 1 (Hollenbach & Tielens 1997, Bernard 
et al. 1993). For this reason, the interiors of (quiescent) 
clouds should be relatively free of UIB emission, and only 
their outer envelopes should significantly contribute to the 
mid-IR background and foreground (e.g. Bernard et al. 
1993, Boulanger et al. 1998). 

As a first step, we neglected the spatial fluctuations of 
the background and foreground, Shack(r) and <5// ore (f), 
and thus assumed hack and If ore to be equal to their mean 
values across the field, < hack > and < If ore > (see 3.5.2 
below for a critical discussion of this assumption). 

For a given gas-to-dust ratio, the H2 column density 
Nh 2 is related to the dust opacity by T\(f) = a\ ■ Nh 2 (r), 
where o\ is the dust extinction cross section^ at the wave- 
length of observation A. We assume a\ to follow the Draine 
& Lee (1984) dust model in the mid-IR as a first approxi- 
mation, but the value is rather uncertain (see section |3.5|) . 
We adopted <T 6 . 75fim = 1.2 x 1CT 23 cm 2 for the LW2 fil- 



ter, 07.75 Mr , 



1.35 x 10 23 cm 2 for the LW6 filter, and 



1.6 x 10 A cm 2 for the LW9 filter, i.e., we calcu- 



4 At the observed mid-IR wavelengths, the scattering compo- 
nent of o\ should be more than ten times smaller than the ab- 
sorption component for a dust composition following Preibisch 
et al. (1993) and can thus be neglected (e.g. Bohren & Huffman 
1983). 



lated cta at the central wavelength of each filter. A more 
accurate value could be obtained by integrating the extinc- 
tion curve over the response of the ISOCAM filters, but as 
we will see in the following, the exact value of a\ has only 
little influence on the derived column density profiles. 

The quantities hack and If or e are unknown, but they 
can be constrained using independent measurements of 
the H2 column density at two different radii. These con- 
straints enable the ISOCAM absorption profile to be con- 
verted into a column density profile. 

3.3.1. Central column density from millimeter continuum 
data 

We used our 1.3 mm continuum maps to estimate the col- 
umn density averaged over the flat central region of each 
core. As millimeter dust continuum emission is optically 
thin, the column density averaged over the flat inner part 
(of radius Rfiat) of the core (see WMA99) may be de- 
rived from the 1.3 mm flux density integrated over the 
same area, S{ l £ nm (cf. AWM96, MAN98): 

N{%Ln = S{ l ^J[rt lat K 1 . 3rnm vm H B 1 .z mrn {T)], (2) 

where n,f lat is the solid angle of the flat inner region, 
Ki.3 mm is the dust opacity per unit mass column density at 
A = 1.3 mm, /i = 2.33 is the mean molecular weight, and 
Bi.z(T) is the Planck function at A = 1.3 mm, for a dust 
temperature T. We assumed a single, representative dust 
temperature Td = 12.5 K for all the cores. Based on recent 
ISOPHOT measurements of L1544, Oph-D, L1689B, and 
other similar starless cores (e.g. Ward- Thompson & Andre 
1999, Lehtinen et al. 1998), we believe that this value of 
Td is likely to be within ±3 K of the true dust temperature 
in most cases. Considering the additional uncertainty on 
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Fig. 4. IRAM 1.3 mm continuum emission maps smoothed to 13" FWHM resolution (contours) superimposed on 
ISOCAM 7 /im absorption image (greyscale). Contour levels are: a) 20, 40, 60 mJy/beam for L1544; b) 30, 50, 70 
mJy/beam for Oph D; c) 20, 40, 60 mJy/beam for L1709A; d) 10, 30, 50 mJy/beam for L1689B; e) 20, 40, 60 
mJy/beam for L310; f) 20, 40 mJy/beam for L328. 
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Fig. 4. (cont'd) Contour levels are: g) 25, 50, 75, 100 mJy/beam for L429; h) 10, 20, 30, 40 mJy/beam for GF5. 



the millimeter dust mass opacity in dense cores, which we 
took equal to «i.3 mm = 0.005 cm 2 g _1 (see, e.g., Prcibisch 
et al. 1993, AWM96, and Kramer et al. 1999), the esti- 
mates of 

N flat 

listed in Table 2 are uncertain by a factor 
of ^, 2 on either side of their nominal values. 

3.3.2. Outer column density from millimeter line data 

The second value of the column density we can use to cal- 
ibrate the ISOCAM column density profile is given by our 
C 18 O(l-0) line observations of the outer parts of the cores. 
Assuming optically thin emission and Local Thermody- 
namical Equilibrium (LTE), we can determine the value 
of the C 18 column density from the following equation: 

N c ^ {cm~ 2 ) = 4.76 x 10 13 * _ 5 °° / TX(v)dv,(3) 

eXp [ rp k ) J — CO 

where T* A is the measured antenna temperature of the 
C 18 O(l-0) line, T& is the gas kinetic temperature in the 
cloud, and dv is in km s _1 (e.g. Rohlfs & Wilson 1996). 
We used Eq. (2) with T k = 10 K to estimate the C 18 col- 
umn density in the outer parts of the core. This method 
would not be reliable towards the center of the core where 
C 18 O(l-0) may be optically thick. Moreover at the centre, 
the possibility of depletion of gas molecules onto grains 
makes the estimation of the column density unreliable 
(e.g. Kramer et al. 1999). We adopted the following re- 
lation between H2 and C 18 column densities (Frerking, 
Langer & Wilson 1982): Nh 2 = 5.53 x 10 6 iV c i8 O + 1.22 x 
10 21 cm~ 2 . The values of the column density Nf^ derived 
in this way, as well as the distance R ou t from the core 
centre at which it was evaluated, are given in Table 0. 
The exact position where -/V^"* was measured is indicated 
by a triangle-symbol on the images of Fig. [j]. The uncer- 
tainty on the H 2 column density derived from molecular 



line data has two major origins: first, the excitation state 
of the C 18 O(l-0) molecules is not accurately known and 
the LTE hypothesis is only an approximation (cf. White 
et al. 1995), and second, the C 18 abundance ratio with 
respect to H2 is also uncertain in our cores. Based on com- 
parisons with more sophisticated calculations performed 
with a Monte-Carlo radiative transfer code (Blinder 1997) 
for a realistic pre-stellar core model, we judge that our 
LTE estimate of the column density in the outer parts of 
each core, N™* , is good to better than a factor of ~ 2. 

3.3.3. Estimates of the background and foreground 
intensities 

Given the values of and N 1 ^ derived from millime- 

ter observations, it is easy to estimate hack and If ore us- 
ing Eq (|J) for f = Rfiat and f = R ou t, and then deduce 
the column density from the mid-IR intensity at any posi- 
tion in the ISOCAM image. The values of hack and Ifore 
that we calculated in this way, using the independent con- 
straints that hack > and If ore > I Z odh ar e listed in 
Table |^. In addition to 'best' values, we give permitted 
ranges for hack an d If ore that were determined using the 
maximum and minimum values of Njg and Njjf ( given 
the factor ~ 2 uncertainty on both these quantities). 

Another approach consists in estimating the values of 
hack and 2/ ore directly from the ISOCAM images. The 
sum hack + Ifore is constrained by the mean mid-IR in- 
tensity < Imir > measured in the map outside the dense 
core (cf. column 7 of Table [lj; this also corresponds to 
the "baseline level" measured on intensity cuts such as 
those of Fig. ||) . Furthermore, one can set straightforward 
lower and upper limits on the foreground intensity since 
Ifore has to be less than the value of the minimum in- 
tensity in the image, I m im and greater than the zodiacal 



13 



Table 2. Core parameters deduced from millimeter observations. 



Name 


Rflat 




Rout 


a rout 


zodi 


f ore 


^back 




(AU) 


2 

cm 


(PC) 


cm -2 


(MJy/sr) 


(MJy/sr) 


(MJy/sr) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


LI 544 


1900 


4.5 xl0 22T 


0.08 


2xl0 2i 


4.2 




0.7i U oi 


Oph D (LW2) 


3400 


5xl0 22 


0.05 


3.2xl0 21 


5.7 




q i+i o 


Oph D (LW3) 










45 


42 1+ 2 ' 6 


q tc+2-0 
y.o_ 2 4 


L1709A 


6400 


4.5xl0 22 


0.18 


3.1xl0 21 


10.5 


n.3i;-j 




L1689B s 


4000 


4xl0 22 


0.16 


2xl0 21 


4.3 




3.811;° 


L1689B w 


4000 


4xl0 22 


0.29 


4.8xl0 21 


4.3 


A 

*"J-0.4 


*- J -0.7 


L310 


8000 


1.7xl0 22 


0.11 


lxlO 21 


7.1 


7 i+0.3 
'• x -1.2 


19.2±^ 


L328 


5000 


3xl0 22 


0.14 


2.4xl0 21 


7.4 


7-4±?. 5 


34.5t°; 9 3 


L429 


2400 


5xl0 22 


0.26 


8xl0 21 


6.8 


6.8+? 2 


9 a + 1.3 
^•°-0.1 


GF5 


5000 


2xl0 22 


0.26 


2xl0 21 


7.2 


9.1±£ 6 


70.8+i 6 5 6 



Notes to Table |: 

(1) Rfiat is the radius of the flat inner part of the core column density profile (see WMA99 and Sect. 
3.4). For L1544, Oph D and L1689B, Rfiat was taken from the millimeter continuum measurements 
of WMA99 and corresponds to the geometrical average of the major and minor axes of the ellipse. 
For L1709A, L310, L328, L429 and GF5, R f i at is the radius of the flat inner part of the ISOCAM 
profile. 

(2) Nfl®* is the column density averaged over the flat inner regio n of t he core of radius Rfiat as 



« ' ±±-2 ~ r- 1 

deduced from our millimeter continuum measurements (see section p. 3.1 ) 



(3) Rout is the r adius at which the C ls O column density was measured in the outer parts of the 
core (see section 3.3.2). 



(4) Njjf is the column density inferred from C 18 O(l-0) observations at Rout (see section 3.3.2). 

(5) I^idi 1 is the zodiacal light emission in the mid-IR (see col. [8] of Table |l|) . 

(6) , (7) I ^ach & Ifore' are the mid-IR background and foreground intensities as deduced from Eq. 
(0) and the values of the centr al and outer column densities, using the additional constraint that 



I fore > I zodi (see section 3.3.3) 



tThe values of N^f and < Nh 2 > 13 » given in Table 2 of WMA99 are too large by ~15 % and 
~ 70 %, respectively. 



emission, I zodi . Hence, I zodl < I fore < I min and, since 

hack + Ifore = < ImIR >, < ImIR > —Imin < hack < 

< I MIR > —Izodi- 

This second method assumes that < Imir > can be 
measured outside the dense core, which is not necessarily 
the case if, e.g., the core extends beyond the limits of the 
ISOCAM map. Moreover, according to Eq. (1), the mini- 
mum intensity in the image is equal to the intensity at the 
maximum of the absorption (/ m ;„ =< hack > -e~ Tmax + < 
Ifore >), which in practice can greatly overestimate the 
value of < If ore > since one typically has T max ^ 1 (see 
below). For these reasons, we preferred to adopt the first 
approach based on millimeter measurements of the central 
and outer column densities (see above). 

In the case of the Oph D dense core, however, it is 
possible to determine accurate values of hack and If or e 
at more than one wavelength using both methods, since 
large-scale (~ 0.8° x 0.8°) ISOCAM images of the p Oph 
main cloud are available (cf. Abergel et al. 1996). This 
provides an important consistency check. At 6.75 /iin, we 
obtain lfj 7 5 e = 5.7±H MJy/sr and I$ a 7 c 5 k = 9.1+1;° MJy/sr 
using millimeter observations (see Table ||), while we 
find 4.6 MJy/sr < lf™ e < 8.5 MJy/sr and 6.3 MJy/sr 

< 7^7| < io.2 MJy/sr using the second method. The 



agreement is thus very good. The background and fore- 
ground intensities can also be estimated at ~ 15 /ini 
(cf. Table |2|). As can be seen in Fig. |=) and |d) be- 
low, they yield a column density profile for Oph D which 
is virtually identical to that derived at 6.75 /*m. Fur- 
thermore, our estimated values {lf ore = 42.1±^ MJy/sr 

and I^ck = 9. 5^2 '4 MJy/sr) are consistent with the idea 
that the mid-IR background is primarily due to the emis- 
sion of UIB carriers while the foreground is dominated 
by the zodiacal light. Indeed, the emission spectrum of 
aromatic particles is such that the ratio of the LW2-filter 
to LW3-filter intensities is ~ 1 (Boulanger et al. 1996), 
which is precisely what we find for the background emis- 
sion in Oph D. By contrast, the 7-15 /im spectrum of the 
foreground emission is rising and compatible with that of 
the zodiacal light (I zo di ~ 5.7 MJy/sr at ~ 7 /jm and 
Izodi ~ 45 MJy/sr at ~ 15 fim). 

The results of these inter-comparisons on Oph D con- 
firm the validity of the millimeter method we used to de- 
termine hack and i/ ore for the other cores. 

We note that for many sources Ifore ~ Izodi in Ta- 
ble ||, suggesting that the zodiacal emission is often the 
main contributor to the foreground. This situation is ex- 
pected if most of the target clouds are illuminated from 
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behind, which is indeed believed to be the case for the 
p Oph cloud, located in front of the Sco OB2 association 
(e.g. de Geus et al. 1989, de Geus 1992). As most of our 
fields were chosen to be nearby and devoid of strong (e.g. 
foreground) infrared sources, the fact that the front sides 
of the clouds appear to contribute negligible foreground 
emission is likely the result of a selection effect. 

3.4- Radial column density profiles 
3.4.1. Derivation of the profiles 

In order to characterize the radial structure of each core 
and reduce the influence of spatial fluctuations in the 
background and foreground, we first derived mean radial 
intensity profiles by averaging the mid-IR emission in the 
ISOCAM images according to the apparent morphology of 
the core. For L1689B, we averaged the intensity over ellip- 
tical annuli of increasing radii separately for a 40°sector in 
the South and a 40°sector in the Western side of the core. 
For L328, OphD and L429, the intensity was averaged 
over circular annuli of increasing radius. Only the South- 
western sector of the profile was considered for Oph D 
and only the Western half for L328. For very elongated, 
filamentary structures (L310, L1709A, L1544, GF5), we 
averaged cuts perpendicular to the main axis of the cores, 
but only kept the Southern half of the absorption inten- 
sity profile (except GF5 for which the North side was 
kept)^|. The sectors over which the profiles were averaged 
(for the elliptically/spherically symmetric cores L1689B, 
L328 and L429) are shown as dotted lines on Fig. [l] as are 
the ranges of cuts which were averaged for the filamentary 
cores (L1544, L1709A, L310, GF5). 

The radial column density profiles N ff2 (f) that were 
derived in this way using our 'best' values for Iback and 
(see above) are displayed in Fig. || as crosses. The 
vertical cross sizes represent statistical error bars on the 
mean radial profile, estimated from the dispersion of the 
values averaged at each radius, divided by the square root 
of the number of values. In addition to these random errors 
at each radius, there is also a global uncertainty on the col- 
umn density profiles resulting from our uncertain knowl- 
edge of Iback and I fore- The dashed curves lying above and 
below the 'best' profiles limit the range of column density 
profiles that are consistent with the mid-IR absorption 
data, given the permitted range of values for Iback and 
I fore in each case (see Table 2). The upper dashed curve 
corresponds to the profile for the smallest value of Iback 
and the greatest value of // ore , while the lower dashed 



5 For Oph D, the Eastern profile does not stretch far enough 
to enable an analysis at large radii (> 200"). For L328, the 
Eastern profile is disturbed by the presence of additional ab- 
sorption features. For L310, L1709A, L1544 and GF5, the pro- 
file we kept was the longest one not affected by the presence 
of supplementary absorption features (or by flatfield noise as 
in L1544). 



curve corresponds to the greatest value of Iback and the 
smallest value of If ore . It is important to stress that there 
is a simple, global relation between each permitted profile, 
NH 2 (r), and the 'best' profile, N H (f), of the type: 



1 



N Ha (r) = — - In 



0~ X 



1 



.4 



D 



where A and B are two parameters depending only on 
N^ 1 and Nf£, and a\ is again the extinction cross- 
section at the mid-IR wavelength A. 

For comparison with the profiles derived from the ISO- 
CAM images, the thin solid curves in Fig. ^| represent 
the column density profiles of finite-sized sphere mod- 
els with uniform density for r < Rfi a t, p oc r~ 2 for 
Rfiat < t < Redge, and surrounded by a medium of uni- 
form column density. The dotted curves plotted in Fig. || 
are Nu 2 oc f _1 power laws (which correspond to the col- 
umn density profile of a singular isothermal sphere at 
10 K). These models as well as other theoretical mod- 
els will be discussed in Sect. 4. (Note that the convolution 
with the ISOCAM psf has no effect on the models beyond 
an angular radius of ~ 6".) 

3.4.2. Results 

It can be seen from Fig. || that all the column density 
profiles present a flattening in their inner regions - though 
this flattening is only marginal in the case of L429. On 
the other hand, the slopes and shapes of the outer parts 
of the profiles appear to vary from core to core and may be 
separated into two broad categories. For four of the cores 
(GF5, L429, L1709A and the Western part of L1689B), 
the profile beyond the flat inner part follows a near power 
law with projected radius f, i.e. Nh 2 oc f , and does not 
show any clear steepening up to the boundary of the im- 
age. By contrast, for the other cores (L328, Oph D, L310, 
L1544 and the Southern side of L1689B), the profile steep- 
ens beyond some radius R m id, until it merges with the am- 
bient molecular cloud. Three cores in our sample present 
"edges" , defined as regions where the column density pro- 
file drops significantly steeper than a Njj 2 oc f -1 power 
law. Globally, we may divide each column density profile 
into four parts: a flat inner region of radius Rff^, followed 
by a near power-law portion of index mi up to a radius 
Rmid ' a re gi°n with a steep slope m 2 from R™£ to R'J® 
that marks the core "edge", and finally a relatively flat 
plateau once the column density has reached the average 
value it has in the ambient molecular cloud. Depending 
on the cases, one or two parts may be absent. 

The values of the core opacities at 7 /im, the column 
densities, as well as the masses, average slopes, and den- 
sity contrast deduced from the ISOCAM profiles of Fig.p 
are given in Table |§|. The peak column density, N^ a , 
was read on the profiles shown in Fig. ||. (In most cases, 
jy-peak ^ jyfjat because central plateau is not strictly 
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Table 3. Core parameters derived from the ISOCAM absorption analysis 



Source 
Name 


(cm" 4 ) 
(1) 




(cm" 2 ) 
(2) 


7/jm 
'max 

(3) 


ra/SO 
"-flat 

(0.01 pc) 

(4) 


TflSO 
(PC) 

(5) 


piso 

(PC) 

(6) 


flat 

(cm" 3 ) 
(7) 


M } l at 

(M ) 
(8) 


M e d ge 

(M ) 
(9) 


mi 
(10) 


Bi2 
(11) 


1/iat 
"out 

(12) 


LI 544 


5.4xl0 22 


6 


OxlO 22 


0.7 


1.4 


0.028 


0.043 


4xlO b 


0.5 


1.9 




_ 44 +u.7 


20 


OphD 


5xl0 22 


6 


6xl0 22 


0.7 


1.6 


0.028 


0.066 


3.0xl0 5 


0.8 


2.0 


-0.8±°i 


9 O+0.7 


80 


L1709A 


4.5xl0 22 


5 


lxlO 22 


0.6 


3.2 


0.18 


>0.32 


1.2xl0 5 


2.4 


>29 


1 n+0.5 
1 - z -0.4 




70 


L1689B 


4.5xl0 22 


5 


8xl0 22 


0.7 


3.2 


0.096 


0.136 


1.4xl0 5 


2.2 


11 


_1 4+0-3 


_4 3+1-8 


50 


(South) 
L1689B 


4.7xl0 22 


6 


OxlO 22 


0.7 


4 


0.28 


>0.40 


1.2xl0 5 


3.5 


>70 


-1 +a4 

l- u -0.8 


50 


(West) 
L310 


3.4xl0 22 


3 


7xl0 22 


0.5 


4 


0.07 


0.1 


9xl0 4 


1.5 


7.7 


i -°-o.i 




>10 


L328 


4.4xl0 22 


5 


8xl0 22 


0.8 


2.5 


0.03 


0.12 


1.8xl0 5 


1.1 


6.4 


_1 7+04 
1 - ' -0.7 


-1 7+ - 4 


>20 


L429 


9.6xl0 22 


1 


2xl0 23 


1.6 


~1.2 


0.35 


>0.5 


6xl0 5 


~0.4 


>60 


-o.9±8:g 


>600 


GF5 


2.7xl0 22 




3xl0 22 


0.4 


2.5 


0.2 


>0.4 


lxlO 5 


0.7 


>30 


-o-9±g:3 




>140 



Notes to Table |: 

(1) Column density averaged over the flat inner part of the core of radius Rff®. 

(2) Peak column density read from the ISOCAM profiles of Fig. [s| 

(3) Opacity at the centre of the core at ~ 7/xm. 

(4) Rff^t is the radius of the flat inner part of the ISOCAM profiles of Fig. [| and is slightly different from Rfi a t (cf. 
Table |) for L1544, Oph D and L1689B. 

(6) Radius of the core edge (when there is one). 

(7) Average density over the flat inner core region, estimated from a piecewise power-law model fitting the profile 
(see text). 

(8) Integrated core mass within a circular area of radius Rff°. 

(9) Integrated core mass within a circular area of radius R^f^,. 

(10) Index of the power-law fitting the column density profile between Rtf^ and R 1 ^^ (Column 4). 

(11) Index of the power-law fitting the column density profile between Rj^^ (Column 4) and R e ^^ e ■ 

(12) Density contrast between the centre and the edge of the core. For Oph D, L1689B and L1709A, the value of 
riout was taken from the large scale 13 CO observations of Loren (1989a), and for L1544, it comes from Butner et 
al. (1995). As no estimate of n out is available for the remai ning cores in the literature, a lower limit for the density 
contrast was estimated from the model described in Section 4.1. 



flat). The average slopes m\ and ni2 of the column den- 
sity profile between Rjf^ and R 1 ^^, and between R^id 
and -Redge> respectively, were derived by fitting segments 
of straight lines to the data. The maximum values of mi 
and ?7i2 were obtained by assuming the maximum possible 
value of Nff^ together with the minimum possible value 
of Njj 1 (taking into account the factor of ~ 2 uncertainty), 
while the minimum values of m\ and m,2 were obtained by 
assuming the minimum possible value of together 
with the maximum possible value of Nf^. 

3.5. Uncertainties, difficulties, limitations 

3.5.1. Dust extinction curve in the mid-IR 

Among the parameters involved in the relationship be- 
tween the mid-IR intensity and the H2 column density is 
the dust extinction coefficient <j\. This paragraph focuses 
on the influence of a\ on the derived column density pro- 
files. The uncertainty on o~\ stems from two main factors. 
First, the ISOCAM filters encompass a certain bandwidth 
over which the value of <j\ varies. In our study, we used for 
<t\ the value for the wavelength at the centre of the filter 



(ie at 6.75 yum for the LW2 filter and 7.75 jum for the LW6 
filter). Considering the rapid variations (steep increase) 
of <7\ towards 8 yxa (eg. Draine & Lee 1984), the values 
we adopted in § |3.3| above could be underestimated by up 
to a factor of ~ 2. Second, the dust properties in dense 
cores are known to depart from those calculated by Draine 
& Lee (1984) for the diffuse interstellar medium (Kriigel 
& Siebenmorgen 1994, see in particular their Fig. 12). In 
environments like dense cores, in which temperatures can 
be as low as ~10 K and densities as high as ~ 10 5 -10 6 
cm -3 , coagulation of particles, fluffiness and the presence 
of ice mantles on the grains can change the value of the 
dust extinction coefficient, depending on the nature of the 
grains (amorphous carbon or coated silicate grains). By 
using the Draine & Lee value, we could thus be underes- 
timating the value of the extinction coefficient by a factor 
of - 2. 

To check how the column density profile was affected 
by the value of a\ , we calculated it for different values of 
<j\ . As the innermost parts of the core are denser and more 
shielded from the interstellar radiation field (and thus 
colder), they are more likely to be affected by grain co- 
agulation and forming of ice mantles, resulting in a higher 
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Fig. 5. Column density profiles of the target pre-stellar cores (crosses), as well as global error bars (dashed curves), 
derived from the ISOCAM images and adjusted so as to match the values obtained in the flat central region with 
1.3 mm continuum data and in the outer region with C 18 O(l-0) line observations (see § 3J3). Theoretical models are 
plotted for comparison (see § 0). 
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Fig. 5. (cont'd) The dotted line is an Nh 2 oc f _1 power-law (corresponding to the profile of a singular isothermal 
sphere at T=10 K), convolved with the ISOCAM psf. The thin solid line corresponds to the column density profile of 
a finite-sized sphere model with uniform density for r < R/iat, P °c r~ 2 for Rfi a t < r < R e dge, and surrounded by a 
medium of uniform column density. 



extinction coefficient than in the low-density outer parts. 
We therefore determined the column density profiles sup- 
posing a linear increase of a\ from 1.35xl0~ 23 cm 2 on the 
edges of the core to 2.7x 10~ 23 cm 2 (which is twice the pre- 
vious value) at the centre of the core. We also varied the 
value of o\ by factors of 2 and 3 to check the influence of 
a high value of o\ on the shapes of the profiles. 



The results of this study show that if o\ varies within 
a factor of ~ 2-3 the column density profiles change only 
marginally. The same holds true if we assume a linear 
gradient in a\ or in Log(a\) from the centre to the edge 
of the cloud (so that in the end a\ varies within a factor of 
2). To summarize, our conclusions regarding the slopes of 
the column density profiles are not sensitive to the existing 
uncertainties on o\. 
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3.5.2. Simulations of non uniform background and 
foreground intensities 

To check the influence of our assumption that the fore- 
ground and background intensities are spatially uniform, 
we carried out simulations of the total infrared intensity 
measured by the camera (and of the column density in- 
ferred) , if the absorbing dense core is embedded in an am- 
bient cloud of uniform column density, itself bounded by 
a mid-IR emitting shell (cf. Fig. 6 - left panel). We con- 
sidered a spherically symmetric core model with uniform 
density for r < Rfi at , p oc r~ 2 for Rf lat < r < R core 
(corresponding to a N oc f _1 column density profile), and 
uniform temperature Ti n ~ 10 K throughout. (This is sim- 
ilar to the simple model compared to the data in Fig. 5 
and discussed in Sect. 4.1 below.) The parent cloud (cor- 
responding to radii R CO re < r < Rdoud) is also assumed 
to be at Ti n ~ 10 K. By contrast, the outer shell (beyond 
Rdoud) is taken to be much warmer at T out ~ 300 K, so 
that its maximum emission is in the mid-IR. Thus, we 
have a thin outer region emitting mid-IR radiation and a 
cold, dense inner cloud/core only absorbing in the mid-IR. 
Such an idealized model of emission by an outer cloud in 
thermal equilibrium at T out ~ 300 K is clearly not realistic 
and only meant to be illustrative of a background medium 
emitting in the mid-IR. Since it is the thin spherical shell 
beyond R c ioud that produces the background and fore- 
ground, Iback and I fore are not uniform in this model but 
present 'limb-brightened' peaks at f = Rdoud (see Fig. || 
- right panel). 

When Rcore is significantly smaller than R c ioud, our 
simulations show that Iback and Ifore vary only slowly 
over the spatial extent of the absorbing dense core (see 
the example of Fig. 6). In such a situation, the column 
density profile derived from the simulated mid-IR inten- 
sity of the model in the same way as for the ISOCAM data 
is almost indistinguishable from the original model profile. 
On the other hand, if R CO re ~ Rcioud so that the peaks in 
hack and If ore occur immediately outside the core, then 
the profile derived under the assumption that Iback and 
Ifore are spatially uniform departs markedly from the true 
column density profile as the temperature discontinuity is 
approached for f ~ Rcioud- In particular, the derived pro- 
file drops sharply near R CO re, he., displays a sharp edge, 
even when the warm outer shell has the same p oc r~ 2 
density profile as the cold inner cloud core (so that there 
is no physical edge in the density structure but only a dis- 
continuity in temperature). This shows that our method 
of converting intensities into column densities enables us 
to detect the existence of a discontinuity, but does not 
necessarily allow us to conclude about its nature. 

In practice, however, it seems reasonable to assume 
that we are in the situation of Fig. 6 and that hack and 
Ifore have nearly uniform values over distance scales of 
several tenths of parsecs. The p Ophiuchi cloud complex 
may be used as an example. In this case, there is con- 



vincing evidence from the IRAS study of Bernard et al. 
(1993) that the mid-IR background emission originates 
from photo-excited PAH-like molecules distributed in a 
sort of shell corresponding to the outer layers of the large- 
scale 13 CO molecular cloud mapped by Loren (1989). This 
geometry is very reminiscent of Fig. 6 with R c ioud on the 
order of a few pc, i.e., much bigger than the rather small 
(Rcore ~ 0.1 pc) cores we are studying here (see, e.g., 
Table ||). This means that the cores are embedded deep 
inside the layers responsible for the mid-IR emission. An- 
other piece of evidence comes from the ISOCAM images 
themselves which do not show strong intensity fluctuations 
on small scales. 



4. Discussion: Comparison with theoretical 
models of core structure 



As pointed out in Sect. B.4, the radial column density 



profiles derived for the pre-stellar cores of our ISOCAM 
absorption sample (see Fig. ||) are generally characterized 
by four different regimes: a) a flat inner region where the 
column density gradient is much flatter than Nh 2 oc r , 
b) a region roughly consistent with Njj 2 oc f -1 , c) an edge 
where the column density gradient steepens quickly with 
distance from core center and gets as steep as Nh 2 oc f~ 2 
or even steeper (see Table ||), until d) the end of the dense 
core is reached and Nh 2 fluctuates about the mean value 
in the parent cloud (i.e., typically N out ~2x 10 21 cm -2 ). 

In the simple case of an infinite spheroidal core with 
a power-law density structure, the radial density gradi- 
ent, p oc r~ p , is related to the column density gradient, 
Nh 2 oc f~ m , by p = to + 1 (e.g. Adams 1991, Yun & 
Clemens 1991). If alternatively, the core radial density 
structure is not represented by a single power law but 
by piecewise power laws, then the preceding relation be- 
tween p and to is no longer strictly exact in the transition 
regions due to convergence effects. However, these effects 
can be quantified in specific cases. For instance, in the 
case of a sphere of uniform density at both small (r < ro) 
and large (r > ri) radii, and having p oc r~ 2 between ro 
and 7"i , the column density profile gets somewhat steeper 
than Nh 2 oc f _1 between ro and n, but remains always 
shallower than Nh 2 oc f~ 14 . For a similar sphere with 



p oc r~ A instead of p oc r between ro and ri, the column 
density profile is always shallower than Njj 2 K r~ 2 ' 2 . We 
conclude that, if the column density profile is observed to 
become steeper than Nh 2 oc f~ 2,2 , then we may safely 
infer that the underlying density profile becomes steeper 
than p oc r~ 3 at some neighboring radius. Note that this 
conclusion remains essentially valid even in the extreme 
case of a finite-sized sphere where the density drops to 
beyond a radius R ou t (cf- Arquilla & Goldsmith 1985, 
Yun & Clemens 1991). In such a case, both the density and 
the column density gradient become effectively infinite at 
Rout- Although the simple relation p = to + 1 would yield 
too large values of p in a narrow range of radii r R ou t , 
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Fig. 6. Schematic representation of the simple model cloud/core system discussed in the text (left panel). The right 
panel shows the simulated mid-IR intensity profile (solid curve) emerging from such a model. The dashed curve 
represents what the intensity profile would be without the absorbing dense core/cloud. The illustrative case shown 
has Rfiat = 0.025 pc, R cor e = 0.15 pc, and R c ioud = 0.3 pc. Note some similarity with the observed profiles of Fig. 3. 



Table 4. Values of reduced \ 2 f° r the various theoretical models 



Source 


Parameterized 
Sphere 


Logotrope 


best fitting SIS 1 


BM95 


CM95 


L1544 


5.9 


8.2 


25.8 


11.3 


5.2 


Oph D 


1.3 


13.4 


23.0 


19.9 


1.3 


L1709A 


3.0 


5.5 


8.8 


3.0 




L1689B s 


1.0 


9.8 


29.4 


1.2 


17.7 


L1689B w 


3.2 


47.2 


241.7 


3.7 


21.1 


L328 


4.3 


30.8 


19.5 


49.7 


6.8 


L429 


19.2 


15.4 


14.3 


5.7 




GF5 


8.0 


4.4 


54.0 


3.3 





Notes: 



The x 2 calculations were restricted to radii smaller than R e d ge (Table |3|) for cores 
with an edge. 

1 The SIS model considered in this Table is the one that fits best the near power-law 
portion of the column density profile. It corresponds to gas temperatures higher than 
10 K. 



the basic conclusion that p becomes infinitely large close 
to R ut would still be correct. 

Characteristics a) and b) above emphasize the fact 
that pre-stellar cores flatten out near their centres and 
are consistent with previous (sub)millimeter continuum 
results (Ward-Thompson et al. 1994; AWM96; WMA99). 
Characteristic c) , implying an outer density profile steeper 
than p oc r~ 3 for three or four of our cores (see Table ||) 
is new and has interesting implications. 

In the following subsections, we discuss the physical 
meaning of the inferred profiles by comparing them with 
the predictions of various theoretical models proposed for 
pre-stellar cores. For each model column density profile, 



NMod(r), we calculated a value of reduced \ 2 according 
to the formula: 



1 



X 



E 



[I(n) - I ba cke-' TMIRNM ^ - Ifo 



-,(4) 



where the sum is over the sampled radii, and is the ob- 
servational rms uncertainty estimated by calculating the 
dispersion of the mid-IR intensities that were averaged at 
radius ri in the image, divided by the square root of the 
number of points at ri (cf. Sect. 3.4). These reduced \ 2 
values which provide an estimate of the "goodness" of the 
various models are listed in Table 0. 
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4--1- Comparison with simple parameterized spherical 
models 

Most of the column density profiles obtained in Sect. |3.4| 
can be reasonably well fitted by finite-sized sphere mod- 
els with piecewise power-law density gradients. The model 
column density profiles shown as thin solid curves in Fig. ^| 
were integrated for truncated spheres of uniform central 
density n c up to a radius Rfiat, and of density decreasing 
as p oc r~ p between Rfiat and an outer truncation radius 
R e dge- Beyond R e d ge , the total ambient column density, 
Nh 2 , is supposed to be uniform. These simple parame- 
terized models provide good approximations to the pro- 
files of polytropic gas spheres (e.g. Chandrasekhar 1967, 
Chieze 1987, McKee & Holliman 1999) and more gen- 
eral equilibrium models for externally-heated, thermally- 
supported self-gravitating spheroids (Falgarone & Puget 
1985, Chieze & Pineau des Forets 1987). The special case 
p = 2 is the most relevant here since the power-law por- 
tion of the observed column density profiles is very close 
to Nh 2 oc f- 1 (see above). In this case, the models mimic 
the column density profiles of pressure-bounded Bonnor- 
Ebert spheres (Bonnor 1956, Ebert 1955). These corre- 
spond to the spherical solutions of the equation of hy- 
drostatic equilibrium for isothermal self-gravitating cores 
confined by external pressure. They have been considered 
as plausible starting points for the gravitational collapse 
of protostellar cores in quiescent clouds (e.g. Shu 1977, 
Foster & Chevalier 1993). 

The parameters n c , Rfiat, and R e d ge , of the 'best-fit' 
models for p — 2 are listed in Table || for each of the cores 
with an edge. These Bonnor-Ebert-like models generally 
yield the lowest values of reduced x 2 among the models 
we have investigated (cf. Table ^). We stress, however, 
that since most of our cores exhibit non-spherical mor- 
phologies they can at best be roughly approximated by 
such models. Furthermore, the density contrasts (~ 10-80 
- see Table ^) that we infer from center to edge for our 
cores are generally larger than the maximum contrast of 
~ 14 for stable Bonnor-Ebert spheres. If thermal pressure 
were the only source of support against gravity, most of 
the cores in our sample should thus be gravitationally un- 
stable. However, since the lifetimes estimated for similar 
starless cores are £ 10 6 yr (e.g. Lee & Myers 1999, Jessop 
& Ward-Thompson 2000), i.e., a factor of £ 10 larger than 
free-fall timescales given the central densities estimated in 
Table it is likely that the cores experience extra sup- 
port from nonthermal pressure forces, especially in their 
outer regions. In this respect, the 'non-isentropic' multi- 
pressure polytropes recently proposed by McKee & Holli- 
man (1999), which include a stabilizing magnetic pressure 
component allowing large density contrasts, offer an in- 
teresting solution (although they do not account for the 
aspherical core morphologies). 



Table 5. Characteristics of the logotropes compared to 
the data. 



Name 


Ps/k 


Pc/Ps 


n c 


Rfiat 






cm" 3 K 




cm -3 


10" 2 pc 


pc 


LI 544 


2.3xlO b 


2.68 


6x10" 


0.3 


0.07 


Oph D 


2.3xl0 6 


2.68 


6xl0 5 


0.3 


0.07 


LI 709 A 


2xl0 5 


20 


2xl0 5 


0.35 


0.4 


L1689B s 


1.4x10" 


4.16 


6xl0 5 


0.3 


0.13 


L310 


2.3xl0 6 


2.68 


5xl0 5 


0.3 


0.07 


L328 


7xl0 5 


4.16 


2.8xl0 5 


0.4 


0.19 


L429 


2.5xl0 5 


20 


5xl0 5 


0.3 


0.33 


GF5 


4.9xl0 5 


4.16 


2xl0 5 


0.5 


0.22 



4-2. Comparison with logotropic spheres 

McLaughlin & Pudritz (1996, 1997) have advocated a 
special non-isothermal equation of state for star-forming 
molecular clouds of the 'logotropic' form P/P c = 1 + 
Aln(p/p c ), where P c and p c denote the central values of 
the pressure and density, respectively, and A is a con- 
stant whose realistic value is A rs 0.2. Assuming this 
equation of state, McLaughlin & Pudritz show that the 
hydrostatic equilibrium solutions for pressure-truncated 
self-gravitating spheres are analogous to Bonnor-Ebert 
isothermal spheres with flat inner regions but with p cx 
r~ l (rather than p oc r~ 2 ) density gradients in their outer 
parts. Stable logotropic equilibria exist for much larger 
density contrasts (up to £ 100) than in the isothermal 
case. 

We compared the column density profiles inferred from 
our ISOCAM data with the profiles of finite-sized lo- 
gotropes at a gas temperature^ of T=10 K, bounded by 
a surface pressure P s . We investigated various values of 
P s and of the center-to-edge pressure contrast P c /P s (see 
Table ||). Once these parameters are fixed, the logotrope 
is entirely defined and the central density, n c , the radius 
of the flat inner region, Rfiat, and the truncation radius 
of the model, R e dge, can all be determined (see Table 1 of 
McLaughlin & Pudritz 1997 and Table §). 

Figure [?] shows (as solid curves) the best-fitting lo- 
gotrope (to which we added a uniform baseline in column 
density) superimposed on the graphs of the experimen- 
tal column density profiles of L1689B (South sector) and 
L328. 

From examination of Fig. ^, it appears that the lo- 
gotropic model cannot match the observed core profiles. 
The disagreement is due to a difference in the overall 
shape/slope of the column density profiles beyond the 
flat inner region: while the observed profiles approach 
Nn 2 (f) oc r -1 , the logotrope is substantially flatter than 
this [with approximately Nu 2 (f) oc ln( 2R * dg * ) for f << 

6 The results of this comparison would remain unchanged 
if we used T=12.5 K, the dust temperature adopted for our 
pre-stellar cores - see Sect. 3.3.1. 
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Fig. 7. Comparison of the logotropic model of McLaughlin & Pudritz (1996) with the column density profiles of 
L1689B (south) and L328. The column density profile of this model is too shallow to fit the observations. 



Redge]- Furthermore, since all realistic logotropes must be 
truncated at R e d ge ~ 0.5 pc to remain stable against ra- 
dial perturbations (see Table 1 of MacLaughlin & Pudritz 
1997), they cannot account for "edgeless" cores such as 
L429, GF5, and L1709A. More quantitatively, the reduced 
X 2 values for the 'best-fit' logotropes are always larger 
than ~ 5, and it is always possible to find a better model 
with a substantially lower value of \ 2 ( see Table ||) . 

4-3. Comparison with self-similar, singular models 

The most extreme versions of the isothermal and lo- 
gotropic models considered above are the singular isother- 
mal sphere (SIS - e.g. Shu 1977) and the singular lo- 
gotrope (e.g. McLaughlin & Pudritz 1997). These are 
hydrostatic self-gravitating spheres with power-law den- 
sity profiles such as p(r) — (a 2 /2irG)r~~ 2 and p(r) = 
(AP c /2nG) 1 / 2 r~ 1 , respectively, where a is the effective 
isothermal sound speed. Although the singular models 
are well into the unstable regime of Bonnor-Ebert and 
logotropic spheres, they present the great advantage of al- 
lowing semi-analytical similarity solutions for the collapse 
(Shu 1977, McLaughlin & Pudritz 1997). Such singular 
models could a priori be representative of the most ad- 
vanced pre-stellar cores, i.e., those on the verge of col- 
lapse and/or protostar formation. In particular, the SIS 
(or slight variants of it) provides the collapse initial con- 
ditions in the widely used 'standard' paradigm of isolated 
star formation (e.g. Shu et al. 1987). 

These singular models, and especially the SIS, are 
much too steep at small radii to match the flat inner region 
characterizing the observed pre-stellar profiles (cf. Fig. B; 
see also Ward-Thompson et al. 1994, 1999, and AWM96). 
Furthermore, as such, these models are scale-free and do 



not account for the presence of edges either. The combina- 
tion of these two qualitative differences with the observed 
profiles imply fairly large values of \ 2 ( see Table ||) . 

The fact that the observed profiles are much flatter 
than the singular models at small radii could imply either 
that the cores of our sample are not sufficiently evolved 
to be described by such models (cf. Li & Shu 1996), or 
that the collapse initial conditions are non-singular. Since 
there is already evidence for significant inward motions 
in (at least some of) our cores (see Sect. 4.4 and Fig. 10 
below), we favor the latter interpretation. 

It is also interesting to note that, even when a large 
part of the observed column density profile is compatible 
with a N H2 oc f" 1 power law (e.g. L429, L1709A), the 
profile of a SIS at T=10 K falls significantly below the ac- 
tual core profile. In order for the SIS profile to match the 
power-law portion of the observed pre-stellar profiles, the 
effective sound speed must be typically twice as large as 
the isothermal sound speed, (a w 0.2 km s _1 for T=10 K). 
This corresponds to an overdensity factor of ~ 4 compared 
to a SIS at 10 K, which again suggests that there is an ad- 
ditional force (e.g. magnetic or turbulent) supporting the 
cores against collapse. In the case of additional turbulent 
support, however, we can exclude the models with ther- 
mal and nonthermal (TNT) support proposed by Myers 
& Fuller (1992), on similar grounds as above, i.e., these 
models are also singular and do not account for the inner 
flattening of the cores. 

4- 4- Comparison with magnetic ambipolar diffusion 
models of cloud cores 

A natural way of accounting for flat inner density 
gradients, high density contrasts, nonthermal support, 
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Fig. 8. Comparison of a magnetically supported core model from BM95 with the south and west column density 
profiles of L1689B. Parameters of the magnetic model are: initial central mass-to-flux ratio /i c o = (M/$) c o = 0.1 (in 
units of the critical value (M/$)^™ t *), initial central density n Cj o = 5.1 x 10 3 cm~ 3 (see also Table ^|) 



and sharp edges is to consider models of magnetically- 
supported cores undergoing ambipolar diffusion (e.g. 
Ciolek & Mouschovias 1994, Basu & Mouschovias 1995 - 
hereafter BM95). According to these models, two phases 
can be distinguished in the evolution of a pre-stellar core. 
The first phase is a long, quasi-static contraction through 
ambipolar diffusion. Along the magnetic field lines, ther- 
mal pressure balances gravitational forces while in the per- 
pendicular directions, the ions are retained by the mag- 
netic field and slow down the process of contraction. Dur- 
ing this magnetically subcritical phase, the central mass- 
to-flux ratio (M/$) cent increases until it reaches the crit- 
ical value (M/$)^* t = (l/2vr)G~ 1 / 2 (e.g. Mouschovias 
1995, BM95). At this point, the inner central core be- 
comes magnetically supercritical and the evolution speeds 
up. The supercritical core collapses dynamically inward, 
while the outer (subcritical) envelope is still efficiently 
supported by the magnetic field and remains essentially 
"held in place". As a result, a (very) steep density pro- 
file or 'edge' develops at the boundary of the supercritical 
core. The magnetic forces induce a disk-like geometry for 
the core, with the polar axis lying along the field lines. 
The volume density of the core is nearly uniform inside 
a roughly spherical central region whose size corresponds 
to the instantaneous Jeans length. Outside this central re- 
gion, the density declines as p(r) oc r~ s , where the power- 
law index s varies typically between ~ 1.5 and ~ 4.5 de- 
pending on radius and on model parameters such as the 
initial mass-to-flux ratio (see, e.g., Fig. 8 of BM95). 

The column density profiles inferred from our data are 
qualitatively consistent with the predictions of these mag- 



netic models, shortly after the formation of the magneti- 
cally supercritical core. 

In particular, we have quantitatively compared our 
(South and West) radial profiles for L1689B to a range 
of theoretical profiles at various evolutionary times cal- 
culated by BM95 (cf. Fig. ||). We find a good agreement 
between the data and the (rotating) core model 7 of BM95 
near time t2, i.e., the time at which the central density has 
increased by a factor of 10 2 (and the central column den- 
sity by a factor 10) compared to its initial value. This is 
close to the start of dynamical collapse in the model. The 
model profiles plotted in Fig. || correspond to a central 
number density of 5.1 x 10 5 cm -3 and a temperature of 
10 K. The Southern and Western profiles of L1689B can 
both be fitted with the same magnetic model if we as- 
sume that L1689B has the thin disk geometry defined by 
BM95 and is inclined by an angle i ~ 65° with respect to 
the line of sight. Thus, the intrinsic column density of the 
core along the magnetic axis, Ng , is a factor cosi ~ 1/2 
smaller than the observed (line-of-sight) column density 
i.e., Ntf 2 ~ 2.3 x 10 22 cm -2 at core centre in this case. 
Furthermore, in the thin disk approximation, the apparent 
minor axis radius of the core (in the north-south direction 
- see Fig. [lji) should also be a factor cosi ~ 1/2 smaller 
than the major axis radius (in the east-west direction). Ac- 
cording to this comparison, L1689B should have just en- 
tered the phase of dynamical contraction following the for- 
mation of the supercritical core. Given that the edge must 
be situated just beyond R e d g e ~ 0.28 pc ~ 5.7 x 10 4 AU 
(cf. Table |^) , we estimate that the radius of the supercriti- 
cal core should be about 0.03 pc ~ 6000 AU (i.e., about 10 
times smaller than R e dge according to model 7 of BM95), 
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which is only slightly larger than the radius of the flat 
inner region (~ 4000 AU). 

An independent argument supporting the claim that 
L1689B is dynamically contracting is provided by the 



L1 689E 



.H 2 CO(2 12 -1„) 




2 4 
Velocity (km/s) 



Fig. 9. H 2 CO(2i2 - In) and N 2 H+(101-012) line profiles 
(solid and dashed histograms, respectively) observed by us 
toward the center of L1689B with the IRAM 30m telescope 

shape of the line profiles we could observe in optically 
thick molecular transitions such as H2CO(2i2 — In) with 
the IRAM 30m telescope (cf. Fig. |^). Toward the center of 
the L1689B core, the H 2 CO(2i 2 - In) line profile is self- 
absorbed and skewed to the blue, while optically thin lines 
(such as the isolated hyperfine component of the N 2 H + (1- 
0) transition) are single-peaked and peak at the dip of the 
H 2 CO self-absorption. This "infall asymmetry" is a classi- 
cal indicator of significant inward motions (e.g. Zhou et al. 
1993, Evans 1999 and references therein). From the shape 
of the H 2 CO(2i 2 — In) line profile shown in Fig. || an in- 
ward speed of ~ 0.07 km s _1 can be inferred for L1689B 
using the simple analytic model of Myers et al. (1996). A 
similar ~ 0.1 km s _1 infall velocity was reported by Tafalla 
et al. (1998) for L1544 over spatial scales of ~ 0.01 pc to 
~ 0.1 pc (see also Williams et al. 1999). 

A variant of the magnetic models considers a non- 
rotating core, in which interactions with the interstellar 
UV radiation field and between neutral and charged grains 
are taken into account (Ciolek & Mouschovias 1994, 1995 
- hereafter CM95). The effect of the UV field is to in- 
crease the degree of ionization in the envelope^, thereby 
improving the magnetic support of the outer enveloppe 
and making ambipolar diffusion less effective there. As a 
result, once a supercritical core has formed, the gap be- 
tween the rapidly infalling central core and the magneti- 
cally supported envelope is enhanced with respect to the 
situation where cosmic rays are the only source of ion- 



7 According to CM95 (see also McKee 1989), UV ionization 
exceeds cosmic ray ionization for A v < 4, which means that 
the UV field can influence cloud support up to A„ ~ 4, while 
it excites UIB carriers only up to A„ ~ 1 (see Sect. 3.S). 



ization. We thus adopted model Byy of CM95 to fit the 
profiles of those of the cores that have an edge (see, e.g., 
Fig.©. 

Using either model 7 of BM95 or model B uv of CM95, 
we were able to find a set of parameters yielding reason- 
ably good fits (reduced % 2 ~ 1-5) for most of our cores 
(see Table ^). We assumed the cores to be thin disks in- 
clined by an angle i with respect to the line of sight and 
applied a scaling factor 7 in column density (and 1/7 in ra- 
dius) to the nominal dimensional models given by BM95 
and CM95, so as to bring the width of these models in 
agreement with that of the cores (while leaving the di- 
mensionless parameters unchanged). 

There are two problems, however, with these ambipo- 
lar diffusion models which involve only a simple static 
magnetic field. First, they produce highly flattened, disk- 
like structures and have difficulties explaining the large- 
scale, filamentary shape characterizing many of the cores 
we observed (e.g. Fig. [j]). A disk- like cloud could appear 
elongated/filamentary if seen nearly edge-on (cf. Ciolek & 
Basu 2000), but viewing most cores edge-on is statistically 
unlikely. The addition of turbulence in the models could 
perhaps alleviate this problem since simulations of turbu- 
lent molecular clouds tend to produce complex filamentary 
structures (e.g. Padoan et al. 1998, Ostriker, Gammie, & 
Stone 1999). Although the molecular lines observed to- 
ward the core centers tend to be rather narrow (see, e.g., 
Fig. § where AV[N 2 H+(101 - 012)] - 0.4 km/s at the 
centre of L1689B; see also Tafalla et al. 1998), indicating 
little turbulence in the core interiors, significant turbulent 
motions are likely to exist in the outer regions, i.e., near 
or beyond R e dge (cf. Goodman et al. 1998). 

Second, and more fundamentally, they require fairly 
large values for the static magnetic field, B re f, in the 
(large-scale) background around the cores. Indeed, it ap- 
pears that ambipolar diffusion models can explain the 
development of sharp edges with s £ 3 only if they are 
highly subcritical initially (see Fig. 8 of BM95). For in- 
stance, model 7 of BM95, which provides a good fit to 
the profiles of L1689B (see Fig. || and Table ^), has an 
initial central mass-to-flux ratio fi C Q = (M/$) c0 = 0.1 (in 
units of the critical value), implying B re f ~ 170 fiG. This 
strong background magnetic field is difficult to reconcile 
with available Zeeman measurements which indicate typi- 
cal field strengths (or upper limits) £ 10-20 [iG on ~ pc 
scales in low-mass star-forming clouds such as Ophiuchus 
and Taurus (e.g. Troland et al. 1996, Crutcher 1999). The 
models taking UV ionization into account (CM95) are 
more satisfactory in this respect, since they can explain 
some of our profiles (L1544, Oph D) while having a back- 
ground magnetic field of only B re f ~ 35 — 50 ^G (cf. 
Table |J), in better agreement with the field values mea- 
sured in dark clouds. Nevertheless, even the CM95 model 
is highly subcritical with an initial central mass-to-flux 
ratio /i c o = 0.256. Such a low value of /x c o must be rare 
in actual cloud cores, since a recent analysis of available 
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Fig. 10. Comparison of the magnetically supported core model Tiuv (at time t2) from CM95 with the column density 
profiles of L1544 and Oph D (see Table || for model parameters). 



magnetic field observations suggests that molecular cloud 
cores are most likely supercritical by a factor of ^ 2 on av- 
erage (Crutcher 1999). Uncertainties in the measurements 
and in the appropriate corrections for projection effects 
still allow critical or marginally subcritical clouds at the 
present time, but values of /z c o significantly smaller than 
0.5 seem unlikely. That cloud cores are not highly subcrit- 
ical, i.e., not strongly magnetized, is further supported by 
the large values of the infall speeds measured in some of 
them such as L1544 (Williams et al. 1999), as well as by 
the slight misalignment between the direction of the mag- 
netic field and the minor axis of the cores observed through 
submillimeter polarimetry (Ward-Thompson et al. 2000). 

It is noteworthy that the ambipolar diffusion model re- 
cently proposed by Ciolek & Basu (2000) for L1544, which 
uses p c0 = 0.8 and B re f rj 12 pG (consistent with Zee- 
man measurements - Crutcher 1999) and matches the cen- 
tral column density distribution of Fig. 5a as well as the 
~ 0.1 km s _1 infall velocities measured by Williams et 
al., does not account for the sharp edge (s Si 4) we infer 
in this core (cf. Table || and Fig. lb of Ciolek & Basu). 

To summarize, although published ambipolar diffusion 
models can potentially explain several of the features we 
observed, it is becoming increasingly clear that they do 
not include all the relevant effects and are not entirely 
satisfactory. In particular, the influence of a turbulent, 
non-static magnetic field on the (outer) structure of pre- 
stellar cores should be investigated more quantitatively in 
the future^. 



. 5. Filamentary structures and cloud fragmentation 



8 When evaluated for the CM95 and Ciolek & Basu (2000) 
models, the cutoff wavelength below which Alfven waves can- 
not propagate in the neutrals (e.g. Eqs. [la] and [2a] of 
Mouschovias 1991) is Xa ~ 0.08-0.14 pc, i.e., comparable to 



As already stated (£4.1), most of the cores observed here 
depart from the spherical morphology adopted in many 
theoretical models (e.g. §4.1 to 4.3). In particular, GF5 



has an elongated, filamentary aspect reminiscent of L977 
and IC5146 discussed by Alves et al. (1998) and Lada, 
Alves & Lada (1999). These authors use extinction mea- 
surements of background starlight in the near-IR to de- 
rive the density structure of their cores. They model L977 
and IC5146 as cylindrically symmetric cores with p oc r~ 2 
density gradients, and point out that their data are not 
consistent with models of isothermal cylinders (which pre- 



dict p 



Ostriker 1964). Lada et al. thus invoke new 



models of cylindrical clouds threaded by helical magnetic 
fields (Fiege & Pudritz 1999a), which yield density pro- 



files falling off as 



-1.8 



to r 



in agreement with their 



extinction results. Our data however (see, e.g., GF5 in 
Fig. and Fig. |Jh) suggest that such elongated cores are 
made up of a succession of discrete spheroidal condensa- 
tions with outer density gradients approaching p oc r~ 2 , 
rather than of a single continuous, cylindrical filament 



with p ' 



' . Our interpretation is supported by the fact 



that the condensations of, e.g., GF 5 and Oph D are seen 
in more than one tracer or wavelength. 

We have derived the (half-power) radii (R ~ 3500- 
12000 AU), mean column densities (< Nh 2 >~ 2- 
6xl0 22 cm~ 2 ) and densities (< nu 2 >~10 5 -10 6 cm~ 3 ) of 
the three condensations of the GF5 filament and of the two 
condensations of the Oph D core and inferred their masses 



Hedge m Table ^. This suggests that hydromagnetic turbulence 
may indeed play a role near or beyond Ridg°e but quickly decays 
on smaller scales. 
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Table 6. Parameters of adopted magnetic models 



Source 


Model 


Inclination angle i 


Scaling factor 7 


n c ,o (cm - ") 


B re f (uG) 


reduced 


L1544 


model B uv of CM95 at t 2 


75 ° 


1 


2.6 xlO 3 


35 


5.2 


OphD 


model B uv of CM95 at t 2 


75 ° 


1 


2.6 xlO 3 


35 


1.3 


L1709A 


model 7 of BM95 at t 2 


70 ° 


1.5 


4.5 xlO 3 


150 


3 


L1689B 


model 7 of BM95 at t 2 


65 ° 


1.7 


5.1 xlO 3 


170 


~ 2 


L328 


model Buv of CM95 at t 2 


° 


1.4 


3.6xl0 3 


50 


6.8 


L429 


model 7 of BM95 at t 3 


65 ° 


1 


3 xlO 3 


100 


5.7 


GF5 


model 7 of BM95 at t 2 


60 ° 


1 


3 xlO 3 


100 


3.3 


Notes: 



(1) t-2 (resp. t'i) is the time at which the central number density has increased by a factor of 10 2 (resp. 10 3 ) 

(2) n c ,o is the initial central density 

(3) B re f is the large-scale magnetic field in the reference state 

(4) u c0 = (M/$)co = 0.1(M/$)^?t for model 7 of BM95 and (5) /x c0 = 0.256(M/$)^£* for model B uv of CM95 

(6) 7 is the scaling factor by which we have to multiply the column density and to divide the radius of the models given by 
BM95 and CM95 so as to adjust the widths of the models to that of our cores. 

(7) The temperature of the cloud is 10 K. 



from our ISOCAM and 1.3 mm continuum measurements. 
The estimated condensation masses, around 0.9 M Q for 
GF5 and 0.3 M for Oph D are consistent with, or slightly 
larger than, the Bonnor-Ebert critical mass (e.g. Bon- 
nor 1956), Mbe — 2ARa 2 s /G (where a s is the isothermal 
sound speed), of 0.7 M Q and 0.2 M©, respectively, for a 
gas temperature of 10 K. This suggests that the conden- 
sations are self-gravitating [M > Mbe)- In the case of 
Oph D for which we have DCO + (2-l) observations, we 
also calculated the virial masses, M vir = 5Ra^ ot /G, of the 
condensations assuming their density is uniform (which is 
consistent with the flattening of the column density profile 
of the cores in their inner part). The measured DCO + (2-l) 
linewidths are ~0.5 km s~ x , yielding M/M V i r ~ 0.5 and 
confirming that these structures are gravitationally bound 
(cf. Pound & Blitz 1993). The characteristics of the GF5 
and Oph D condensations are reminiscent of those found 
for the pre-stellar fragments of the p Oph protocluster 
(MAN98). The average spacing between two condensa- 
tions is 20000 AU for GF5 and 9000 AU for Oph D, which 
is of the order of the Jeans diameter Aj ~ 12000 AU in 
the parent cloud core, for a gas temperature of 10 K and 
an average density of 10 5 cm -3 . Furthermore, the ratio 
of condensation spacing to filament diameter is ~ 2, con- 
sistent with models where gravity plays a major role in 
fragmenting the cores/filaments (Larson 1985, Nakamura 
et al. 1993, Fiege & Pudritz 1999b). 

The existence of these compact condensations also em- 
phasizes the fact that even isolated prestellar cores (such 
as most cores in our sample) are fragmented and are likely 
to form more than one star, perhaps even small groups of 
stars (cf. Tafalla et al. 1999). 



suggests that the initial conditions for fast protostellar 
collapse are often neither singular nor scale-free. In many 
cases, the collapse initial conditions probably differ sub- 
stantially from the idealized singular isothermal sphere 
(SIS) adopted by the standard theoretical scenario of Shu 
et al. (1987) which leads to a constant mass accretion rate 
during the protostellar phase. It is likely that protostar 
evolution will depart, at least initially, from the predic- 
tions of the standard paradigm, even if this model provides 
a good, first-order approximation of self-initiated, iso- 
lated star formation. When starting from Bonnor-Ebert- 
like initial conditions, all (magneto-)hydrodynamic col- 
lapse models predict a decline of the accretion rate dur- 
ing the protostellar phase (e.g. Foster & Chevalier 1993, 
Tomisaka 1996, Henriksen et al. 1997, Basu 1997, Safier, 
McKee, & Stahler 1997, Ciolek & Konigl 1998, Li 1999). 
Interestingly enough, observational evidence for a declin- 
ing protostellar accretion rate has been found in the form 
of a marked decrease of outflow powei|] with time from 
young (Class 0) protostars to evolved (Class I) protostars 
(Bontemps et al. 1996). 

Another related implication of the sharp outer edges 
we observed is that most pre-stellar cores appear to be 
effectively decoupled from their large-scale parent clouds, 
providing a finite reservoir of mass for subsequent star for- 
mation. This supports the view according to which there 
is a mass scale in the star formation process (e.g. Larson 
1999) and the stellar initial mass function is at least partly 
determined by fragmentation at the pre-collapse stage (see 
also MAN98). 



4-6. Implications for the earliest stages of collapse 

The presence of flat inner radial density gradients and 
sharp outer edges in the pre-stellar cores of our sample 



Since protostellar outflows are believed to be directly pow- 
ered by accretion (e.g. Konigl & Pudritz 2000), the outflow 
momentum rate may be used as a surrogate tracer of the ac- 
cretion rate. 
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5. Summary and conclusions 

We imaged a sample of 24 pre-stellar cloud cores at 7 pm 
with the mid-IR camera ISOCAM aboard the ISO satel- 
lite. Our main results and conclusions are the following: 

1. As many as 23 of the 24 fields of our programme show 
absorption-like features in the mid-IR. By observing 
and detecting 8 of these in emission at 1.3 mm with 
the IRAM 30 m telescope, we show that most of the 
mid-IR dark features trace dense, cold core material 
seen in absorption against a diffuse background emis- 
sion arising from the rear side of the parent molecular 
cloud. 

2. Based on a simple analysis of the mid-IR intensity pro- 
files, and using independent constraints on the central 
and outer column densities from millimeter continuum 
and line data, we are able to derive the mean radial 
column density profile of each core detected in absorp- 
tion with ISOCAM. 

3. We confirm that pre-stellar cores have density and col- 
umn density profiles that flatten out near their cen- 
ters, a conclusion already reached by Ward-Thompson 
et al. (1994, 1999) and Andre et al. (1996) on the 
basis of (sub) millimeter emission maps. At radii less 
than R f i at ~ 4000-8000 AU (i.e. ~ 0.02-0.04 pc) 
the derived radial profiles are significantly flatter than 
p oc r~ 2 and Nh 2 °c r~ x (i.e., as flat as Nh 2 oc f-° 15 - 
f~ - 2 ), respectively. 

4. We find that, beyond R e dge ~ 0.05-0.3 pc, at least 
three (and possibly four) of the cores in our sample 
are characterized by sharp edges corresponding to ra- 
dial column density profiles that drop steeper than 
Nh 2 oc f~ 2 with projected radius. Depending on the 
details of deprojection effects, this suggests that the 
density gradient becomes itself steeper than p oc r~ 3 
near R e dge- An important implication is that the mass 
reservoir available for star formation in these cores is fi- 
nite, supporting the idea that stellar masses are partly 
determined at the pre-stellar stage. 

5. Between Rfi a t and R e dge, the core radial profiles ap- 
proach a near power-law regime consistent with Nh 2 oc 
f _1 and p oc r~ 2 . In particular, we show that models 
predicting much flatter power-law radial density pro- 
files (e.g. p oc r _1 ) such as logotropic non-isothermal 
spheres (e.g. McLaughlin & Pudritz 1996) cannot ac- 
count for the structure of our low-mass pre-stellar 
cores. 

6. Points (3) and (4) above suggest that the initial condi- 
tions for fast protostellar collapse are neither singular 
nor scale-free. Semi-analytical collapse models based 
on self-similar initial conditions such as the singular 
isothermal sphere or the singular logotrope can thus 
provide only a first-order, approximate description of 
the protostellar phase. 

7. We have compared our results with various theoretical 
models of core structure, and have found that ambipo- 



lar diffusion models of magnetically-supported cores 
(e.g. Ciolek & Mouschovias 1995, Basu & Mouschovias 
1995) are most promising although they require a 
(probably unrealistically) strong background magnetic 
field (~ 30-100 pG). We hypothesize that more elabo- 
rate versions of these models, incorporating the effects 
of a non-static, turbulent magnetic field in the outer 
parts of the cores would be more satisfactory and could 
also account for the filamentary shape often seen on 
large ( ^ 0.25 pc) scales. 
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